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Abstract 

We study the possibility that the 7-ray loud active galactic nuclei (AGN) are the sources of 
ultra-high energy cosmic rays (UHECR), through the correlation analysis of their locations and 
the arrival directions of UHECR. We use the 7-ray loud AGN with d < 100 Mpc from the second 
Fermi Large Area Telescope AGN catalog and the UHECR data with E > 55EeV observed by 
Pierre Auger Observatory. The distribution of arrival directions expected from the 7-ray loud AGN 
is compared with that of the observed UHECR using the correlational angular distance distribution 
and the Kolmogorov-Smirnov test. We conclude that the hypothesis that the 7-ray loud AGN are 
the dominant sources of UHECR is disfavored unless there is a large smearing effect due to the 
intergalactic magnetic fields. 
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I. INTRODUCTION 



The origin of ultra-high energy cosmic rays (UHECR), whose energies are above 1 EeV(= 
10^^ eV), has been searched for many years; however, it is still in vague. In the search for the 
origin of UHECR, the Greisen-Zatsepin-Kuzmin (GZK) suppression jl, 2| plays an important 
role. This suppression tells us that the sources of UHECR with energies above the GZK 
cutoff, Eqzk ~ 40 EeV, should be located within the GZK radius, tgzk ~ 100 Mpc, because 
the UHECR coming from beyond the GZK radius cannot reach us by loosing energies as 
consequences of the interactions between cosmic microwave background photons. The recent 
observations support the GZK suppression, thus we can focus on the source candidates 
lying within ~ 100 Mpc. 

For the possible sources of UHECR, several kinds of astrophysical objects have been 
proposed, such as 7-ray bursts, radio galaxies, and active galactic nuclei (AGN) which 
are known to be able to accelerate UHECR enough. To verify that these astrophysical objects 
could be the sources of UHECR, it is worthwhile to compare the arrival directions of UHECR 
and the positions of source candidates u sing statistical tests. Although many statistical 



studies for correlation have 
In our previous work 21 



3een done 



14 



22l |. the origin of UHECR is not confirmed yet. 



22|, we examined the AGN model, where UHECR with energies 



above a certain energy cutoff are generated from the AGN lying within a certain distance 
cut. We used AGN listed in Veron-Cetty and Veron (VCV) catalog [231, bj] and the UHECR 
data observed by Pierre Auger Observatory (PAO) 19|, l20||. By statistical test methods, we 
concluded that the whole AGN listed in VCV catalog cannot be the true sources of UHECR 
and pointed out that a certain subset of listed AGN could be the true sources of UHECR. 
Some subsets of AGN, which have the counterpart in X-ray or 7-ray bands have been tested 
for the possibility that they are responsible for UHECR 
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In the case of the 



correlation with AGN detected in hard X-ray band, it is found that the fractional excess of 
pairs relative to the isotropic expectation 20|. In the case of AGN detected in 7-ray band, a 
marginal correlation is found when we consider the small angular separations only, but the 
strong correlation is found when we consider the angular separations up to 20° [27]. 

This paper focuses on the statistical tests for the correlation of UHECR with certain 
subsets of AGN, especially 7-ray loud AGN because 7-ray loud AGN have sufficient power to 
accelerate UHECR. In Section |Tll the source models of UHECR which assume that UHECR 
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with E > Ec come from AGN with d < dc which are emitting strong 7-rays is presented. 
Based on the 7-ray loud AGN model we construct, we create the mock arrival directions of 
UHECR by Monte-Carlo simulation. In Section UTTl we describe the UHECR data and the 7- 



ray loud AGN data we use in our analysis. We use 2010 PAO data [20|] for observed UHECR 
data and the second catalog of AGN detected by the Fermi Large Area Telescope (LAT) 29 1 
for the 7-ray loud AGN data. To test the correlation between the observed UHECR and 
the 7-ray loud AGN, the methods which can compare the distribution of observed UHECR 
arrival direction and that of the mock UHECR arrival direction expected from the source 
model are needed to be established. The brief descriptions of our test methods are provided 
in Section HVl The results of statistical tests are given in Section |V] and the discussion and 
the conclusion follow in Section IVll 



II. SOURCE MODEL OF UHECR 



Several astrophysical objects are know to be able to accelerate CR up to ultra-high 
energy. Among them, AGN are the most popular objects since the strong correlation was 
claimed by PAO [l^. However, its updated analysis results and other studies exclude the 
hypothesis that the whole set of AGN is responsible for the UHECR 18|, I20l422l|. In our 
previous work 22! , using PAO UHECR having energies above 55EeV {20] and AGN within 
100 Mpc listed in the 13th edition of VCV catalog 2J], we concluded that we can reject 
the hypothesis that the whole AGN within 100 Mpc are the real sources of UHECR. Also, 
we tested the possibility that the subset of AGN is responsible for UHECR; we took AGN 
within arbitrary distance band as source candidates and found a good correlation for AGN 
within 60 — 80 Mpc. However, we do not have a reasonable physical explanation for distance 
grouping. This motivates us to try the subclass of AGN with proper physical properties 
appropriate for UHECR acceleration. 

In this paper, we study the hypotheses that AGN emitting strong 7-ray are the sources of 
UHECR, based on the theoretical study in Ref . 30| . Dermer et al. calculated the emissivity 
of non-thermal radiation from AGN using the first Fermi LAT AGN catalog (ILAC) 31 1 
to confirm that they have sufficient power to accelerate UHECR. In the Fermi acceleration 
mechanism using colliding shell model, they found that some of AGN listed in ILAC have 
enough power to accelerate UHECR. (See the Figure 3. in |30|.) Therefore, we set up the 7- 
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ray loud AGN model 
in our previous work 



or the UHECR source in the same way as our AGN model introduced 



22|. 



When UHECR propagate through the universe, UHECR undergo deflection of its tra- 
jectory by intergalactic magnetic fields. These phenomena are embedded in the simulation 
for the mock UHECR expected from the AGN model to compare with observed UHECR by 
introducing the smearing angle parameter (Og) and by restricting the distance of the source 
(dc) and the energy of observed UHECR (Ec) following the GZK suppression. 

We study two versions of 7-ray loud AGN models in this paper. The first one assumes 
that UHECR with energies E > come from AGN which are emitting strong 7-rays 
with distance d < d^, and the second one assumes that among the 7-ray loud AGN those 
having TeV or very high energies 7-ray emission are responsible for the UHECR. The same 
constraints of UHECR energy and AGN distance are applied to these models. We call the 
first one the 7-ray loud AGN model, and the second one the TeV 7-ray AGN model. 

The UHECR flux in the simulation for these two models are described below. The 
expected UHECR flux at a given arrival direction r is composed of the 7-ray loud AGN 
contribution and the isotropic background contribution 

i^(r) = i^AGN(r) + i^iso, (1) 

where -FAGN(r) is the distribution of all AGN within the distance cut and Figo the contribu- 
tion of isotropic background from the outside of distance cut d^. That is, a certain fraction 
of UHECR is coming from AGN and the remaining fraction of them is originated from the 
isotropic background. We introduce the AGN fraction parameter as 

r _ Fagn , ^ 

~ ^ Tp^' 

-r AGN + -riSO 

where -Fagn = (4vr)~^ / E\Qj<f{r)dQ is the average AGN-contributed flux. 

In the next step, we consider two approaches for -FagnI?) because the relation between 
UHECR flux and AGN property is not established yet. The UHECR flux from all AGN can 
be written as 

FAGN(i-)oc Yl 4^-e^P[-(%(r)/M'] , (3) 

ieAGN 3 

where Lj is the UHECR luminosity, di is the distance, Oj{r) = cos^^(r ■ f^) is the angle 
between the direction r and the j-th AGN, Ogj is the smearing angle of the j-th AGN. 



The first approach assumes that all AGN have the same UHECR luminosity, Lj = L, and 
the same smearing angle, Ogj = Og. The second approach assumes that the UHECR flux 
contributed by AGN is proportional to the 7-ray flux of AGN. 



i^AGN(f)oc F^,,-exp[-(%(f)/e,,)'] 

jeAGN 



(4) 



where F^j is the photon flux of AGN detected by Fermi LAT in the 1 — 100 GeV energy 
band. Although the normalization is needed for the accurate expression of Eq. ([3]) and ([1]), 
we neglect it because we are not concerned with the total flux of UHECR in the test. 

We have two free parameters in our model for the simulation, the smea ring angle Og and 

;he AGN fraction For the flducial values of Og and /a, we take Og = 6° |32| and = 0.7 

33|. 

In the last step for realizing the mock UHECR in the simulation, we need to consider the 
exposure function reflecting the efficiency of the detector. The geometric efficiency of the 
detector depends on the location of the experimental site and the zenith angle cut. Then, 
the exposure function h{6) is given by [34 1 

1 



hi6) 



TX 



[sin a^n cos A cos 5 + sin A sin 5] 



(5) 



where A is the latitude of the detector array, 6^ is the zenith angle cut, and 



ar. 



for e > 1, 

for ^ < -1, with ^ 



cos 9m — sin A sin 5 

cos A cos 5 



cos ^ ^, otherwise 



The latitude of the PAO site is A = —35.20° and the zenith angle cut of the released data is 



Or, 



60° 



III. DESCRIPTION OF THE DATA 



We get the information on the 7-ray loud AGN from the second catalog of Fermi LAT 
AGN (2LAC) published in 2011 y. The 2LAC includes the AGN information collected 
by the Fermi LAT for two years. It contains 1017 7-ray sources located at high galactic 
latitude > 10°) as the low galactic latitude region is masked by the galactic plane. That 
region is excluded because the low galactic latitude region is too noisy to be investigated due 
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Name 


I 


b 


z 


I 


TeV flag 


class 


Centaurus A 


309.52 


19.42 


0.0008 


3.03 X 10-9 


Y 


Radio Galaxy 


NGC 0253 


97.37 


-87.96 


0.0010 


6.2 X 10-1° 


N 


Starburst Galaxy 


M82 


141.41 


40.57 


0.0012 


1.02 X 10-9 


N 


Starburst Galaxy 


M87 


283.78 


74.49 


0.0036 


1.73 X 10-9 


Y 


Radio Galaxy 


NGC 1068 


172.10 


-51.93 


0.0042 


5.1 X 10-1° 


N 


Starburst Galaxy 


Fornax A 


240.16 


-56.69 


0.0050 


5.3 X 10-1° 


N 


Radio Galaxy 


NGC 6814 


29.35 


-16.01 


0.0052 


6.8 X 10-1° 


N 


Unidentified 


NGC 1275 


150.58 


-13.26 


0.018 


1.88 X 10-^ 


Y 


Radio Galaxy 



TABLE I: The 8 clean 7-ray loud AGN within 100 Mpc. /: galactic longitude (degrees), b: galactic 
latitude (degrees), z: redshift, F^: photon flux (photon/cm^ /s), TeV flag: TeV AGN, class: optical 
class 



to diffuse radio emission, interloping galactic point sources, and heavy optical extinction. 
There are 886 AGN samples, which are called clean AGN, categorized by the condition that 
the sole AGN is associated with the 7-ray source and has the association probability P is 
larger than 0.8. 

In the 7-ray loud AGN model, we pick up the distance cut dc = 100 Mpc corresponding to 
the redshift z ~ 0.024. (We use h = 0.70, fi^ = 0.27, and Ha = 0.73 to convert the redshift 
to the distance.) Then, only 8 AGN among the clean AGN are picked up as UHECR source 
candidates. For the TeV AGN model, we use the list of TeV AGN detected by the Fermi 
LAT. There are 34 TeV AGN among the clean AGN and only 3 TeV AGN are used for the 
source candidate after the distance cut. (See the Table 9 in (29j.) In the Table [U we list 
8 source candidate AGN within 100 Mpc, which are detected in 7-ray range, and the TeV 
AGN are marked using TeV flag, Y. 

We use PAO 2010 data [20] for the observed UHECR data, which were collected by the 
surface detector from 2004-01-01 to 2009-12-31. It includes 69 events in the declination band 
S = — 90°-24.8° in the equatorial coordinates and having energies above 55EeV. Among 
them, we use only 57 PAO UHECR to avoid the galactic plane region \b\ < 10° as mentioned 
above. Fig. [1] shows the distributions of the arrival directions of PAO data and 7-ray loud 
AGN detected by Fermi LAT in the galactic coordinates using the Hammer projection. 
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FIG. 1: Distributions of 8 Fermi LAT AGN within 100 Mpc (blue and red squares) and 57 PAO 
UHECR (black bullets) in galactic coordinate using the Hammer projection. The red squares 
represent 7-ray emitting AGN in TeV band. The magenta line means a boundary of PAO field of 
view and the cyan lines represent the border of the low latitude region (j6j < 10°). 

IV. STATISTICAL TEST METHOD 



To compare the arrival direction distribution expected from the model and that of the 
observed data, we need the method by which we can represent the characteristic of the 



distribution of arrival direction and apply t 



comparison methods in the previous paper 



le statistical test easily. We proposed some 
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22| . In this analysis, we take the correlational 



angular distance distribution (CADD) method which is most appropriate for the test of the 
correlation between point sources and UHECR. CADD is the distribution of the angular 
distances between all pairs of the point source and UHECR arrival directions: 



CADD : \^9ij' = arccos ( ■ r^- 



M} 



(6) 



where f , are the UHECR arrival directions, r'j are the point source directions, and and M 
are their total numbers, respectively. Now, we get two CADDs to compare: CADDq from 
the observed UHECR and the 7-ray loud AGN, and CADDm from the mock UHECR of the 
model under consideration and the same 7-ray loud AGN. The total number of the data in 
CADD is A^cADD = NM, which means that the number of data in CADD are larger than 
the sampling number N. 
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By comparing CADDq and CADDm, we can test whether our models are suitable to 
describe the observation. There are several statistical test methods which can prove that 
two distributions are different or not. One of the widely used statistical test methods is 
Kolmogorov-Smirnov (KS) test. It uses the KS statistic, which is the maximum absolute 
difference (-Dks) between two cumulative probability distributions (CPD), CPD of the ob- 
served CADDo, So{x), and that of the theoretically expected CADDm, Sm{x), 

Dks = max \So{x) - Sm{x) \ . (7) 

Once we calculate the KS statistic, we can get the probability that two different distributions 
come from the same population through the Monte-Carlo simulation. To get CADDm accu- 
rately, we generates 10^ mock UHECR data. To obtain the probability distribution of the 
KS statistic -Dks, we generate 10^ -Dks for a given model. Thus our probability estimation 
is reliable up to roughly 10~^. 

V. RESULTS 

In this work, we test 4 models for the UHECR sources: 1) the 7-ray loud AGN model 
with UHECR flux proportional to the inverse square of the distance, 2) the 7-ray loud AGN 
model with UHECR flux proportional to the 7-ray flux of AGN, 3) the TeV AGN model 
with UHECR flux proportional to the inverse square of the distance, and 4) the TeV AGN 
model with UHECR flux proportional to the 7-ray flux of AGN. From now on, we call them 
7-d model, 7-f model, T-d model, and T-f model, respectively. 

Fig. [2] shows the distributions of the mock UHECR of 4 models with the smearing angle 
6s = 6° and the AGN fraction Ja = 0.7. The two left panels are for the 7-ray loud AGN 
model, i.e., 7-d model (upper panel) and 7-f model (lower panel), and the two right panels 
are for the TeV AGN model, i.e., T-d model (upper panel) and T-f model (lower panel). 
The blue squares mark the locations of 7-ray loud AGN and the red dots represent mock 
UHECR generated from the source models. The mock UHECR concentrated near the AGN 
are generated from the source AGN and the uniformly distributed mock UHECR com from 
the isotropic background. 

In Fig. [21 we can see the distinguished features depending on the source models we 
assumed. There are 6 7-ray loud AGN in the fleld of view of PAO (7-d model and 7-f 
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FIG. 2: Distributions of the mock UHECR with smearing angle 9s = Q° and AGN fraction = 1 
(upper panels) and that of the mock UHECR with AGN fraction /a = 0.7 (lower panels). Left 
panels are for the 7-ray loud AGN model and right panels are for the TeV AGN model. The red 
dots are the mock UHECR generated by each AGN model and other marks are same with Fig. [TJ 

model), and there are 2 TeV AGN among them (T-d model and T-f model). Also, we can 
see the difference in the mock UHECR distributions due to the difference in UHECR flux 
modeling. For the 7-d model, Cen A and NGC 0253 are the dominant sources because 
they are close to us. In contrast, for the 7-f model, 6 AGN contribute rather equally to the 
generation of mock UHECR. For the T-d model, Cen A is a dominant source, but Cen A 
and M87 share the proportion to generate the mock UHECR in the T-f model. 

In Fig. [3l the probability of each model as a function of the AGN fraction from to 1 
and the smearing angle from 0° to 180° is given. The black gradient color represents the 
probability that the arrival direction distribution of the PAO UHECR comes from the given 
AGN source model. The red, green, and blue contours represent la, 2a, and 3a probability 
lines. 

First, let us compare the results of the 7-ray loud AGN models and TeV AGN models 
with UHECR flux proportional to the inverse square of the distance. Because the set of 
TeV AGN is a subset of 7-ray loud AGN, the principal difference between two models is the 
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AGN fraction (f) AGN fraction (f) 

FIG. 3: Probability dependencies on the AGN fraction {/a) and the smearing angle Og- The left 
panel is for 7-ray loud AGN model and the right panel is for TeV AGN model. The black gradient 
color means the probability and the solid lines represent the contour plot for la (red), 2a (green), 
and 3a (blue). 

number of source AGN. It brings a significant difference in the results, which is shown in the 
upper panels in Fig. E] (7-d model and T-d model). As a rule of thumb, we choose the 3a 
probability as a criterion for ruling out the model. The critical value of the AGN fraction 
is fA,c ~ 0.4 and the critical value of the smearing angle is 6s.c ~ 25° for the 7-ray loud 
AGN model. Decreasing the AGN fraction from the critical value or increasing the smearing 
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angle from the critical value increases the probability. In comparison, for the TeV model, the 
critical value of the AGN fraction is /a,c ~ 0.2 and the critical value of the smearing angle 
is ds^c ~ 110°. We can deduce that it is hard to describe the observed UHECR distribution 
from T-d model. 

Next, let us compare the results of two models with the same source AGN but with 
different UHECR flux models: the model which assumes that UHECR flux is proportional 
to the inverse square of the distance of source AGN and the model which assumes that 
UHECR flux is proportional to the 7-ray flux of source AGN. When we compare 7-d model 
to 7-f model, the critical values of the AGN fraction and the smearing angle are /a,c ~ 0.4 
and ds^c ~ 25° for the 7-d model, and the critical values of the AGN fraction and the 
smearing angle are fA,c ~ 0.7 and Og^c ~ 20° for the 7-f model. We can exclude the AGN 
models in the cases having the parameters within the critical values (blue line). There is 
no crucial distinction between the different UHECR flux models for the 7-ray loud AGN 
models. 

However, when we compare T-d model to T-f model, the different UHECR flux models 
result in the meaningful distinction. In the T-d model, Cen A is a dominant source to 
generate the mock UHECR, thus they are clustered around Cen A. Nevertheless, there is no 
significant difference among Cen A, M87, and NGC 1275 to produce the mock UHECR in 
the T-f model, thus the mock UHECR are clustered not only around Cen A but also around 
M87 and NGC 1275. This is the reason why the probability is increased dramatically as the 
smearing angle increases in the T-f model. The proportion of the mock UHECR produced 
by NGC 1275 which is located outside of the field of view increases as the smearing angle 
increases. Therefore, the clustered feature of the T-d model is quite different from the T-f 
model and the small number of source cause the clear discrepancy. 

In short, we find that the source models assuming 7-ray loud AGN are responsible for 
UHECR (7-d model and 7-f model) are more plausible than the source models assuming 
AGN having higher energy are responsible for UHECR (T-d model and T-f model). Also, 
which flux model is appropriate for describing the UHECR flux is not conclusive yet. This 
seems worthwhile to continue to study. At this stage, we can state the critical values that 
the null hypotheses are rejected. The critical regions are inside the 3cr contours. For the 7-d 
model, the critical values of the AGN fraction and the smearing angle are /a,c ~ 0.4 and 
Og^c ~ 25°, and for the 7-f model, the critical values of the AGN fraction and the smearing 
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angle are /a,c ~ 0.7 and 9s,c ~ 20°. The critical values are fA,c ~ 0.2 and 6s,c ~ 110° in the 
case of the T-d model and the critical values are /a,c ~ 0.6 and Og^c ~ 30° in the case of the 
T-f model. That is, the 7-ray loud AGN dominance models with small smearing angle are 
excluded. 



VI. DISCUSSION AND CONCLUSION 

The purpose of this work was to test the possibility of a subclass of AGN which emit 
strong 7-rays as the source of UHECR. We cannot confirm that the source model using 7-ray 
loud AGN is better than the source model using whole AGN in explaining the distribution 
of arrival direction of UHECR. Compared to the results of our previous work [22|, which 
assumes that whole AGN in the VCV catalog are the source of UHECR, the contour of 
la is changed and the overall probability distribution for each hypothesis is higher than 
the previous results. (See the probability plot for the 7-d model and the left panel of Fig. 



7 in 



22|.) However, we cannot tell 7-ray loud AGN model describes the observation well 



definitely because the higher probabilities mean that the simulated distribution by models 
;he observed distribution only. 



are consistent with 

Nemmen et al. 27| investigate the correlation of UHECR with 7-ray loud AGN using 
the first 27 PAO data and ILAC. Since the data set and the constraints used in our work 
and those used in their work are different, we need to compare the results carefully. To test 
cross- correlation between PAO and ILAC, they count the cumulative number of correlated 
UHECR as a function of angular distance from the ILAC sources located within 200 Mpc. 
They find that the angular distance ip ^ 17° which minimizes the probability that the 
observed distribution would be the isotropic one. Although we pointed out that this method 
is not suitable to test the cross-correlation directly because it focuses on the deviation from 



2l|, the result is consistent 



isotropy rather than the direct correlation in the previous paper 
with our results that 7-ray loud AGN dominance models having small smearing angle could 
be rejected. 

The study of arrival direction of UHECR and its source could be one way of estimating 
the magnitude of the intervening magnetic fields or identifying the primary particle. If it is 
true that the 7-ray loud AGN could be the source of UHECR, the large smearing angle may 
imply the large deflections by intervening magnetic field. Interestingly, the larger deflection 
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angles than traditional pred.cfon ones are estimated by Ryu et al. M. also. These may 

support the measurements of the shower maximum by PAO [37] — the primary particle is 
presumed to be heavy nuclei. In addition, the results by Dermer et al. 30| show that heavy 
nuclei are more likely to be accelerated to ultra-high energy in AGN. Taken together, it is 
possible that the primary particle would be heavy nuclei. 

In summary, we tested the possibility that the 7-ray emitting AGN are the sources of 
UHECR. We took two sub-classes of AGN. One is the set of AGN emitting strong 7-rays 
observed by Fermi LAT and the other one is the subset of the former having 7-rays in TeV 
band. For the UHECR flux, we considered two possibilities: UHECR flux is proportional 
to the inverse square of the distance times the luminosity of its source and UHECR flux 
is proportional to the photon flux of AGN. We used UHECR with E > 55EeV and based 
on GZK suppression, we restricted AGN to be within dc = 100 Mpc. Also, we introduced 
two free parameters in the simulation to produce the mock UHECR, AGN fraction and 
smearing angle 6s so that we could manipulated the contribution of the source candidate 
and the deflection by intervening magnetic fields. To compare the distributions of observed 
UHECR to that of the mock UHECR, we adopted CADD which is suitable to test cross- 
correlation between UHECR and sources. The probabilities for 4 models calculated by 
KS test using CADD tells us that TeV AGN model with UHECR flux proportional to the 
inverse square of the distance is not appropriate to depict the observation. Rather, TeV AGN 
model with UHECR flux proportional to the photon flux is more plausible to describe the 
observation. In the cases of 7-ray loud AGN models, the effects of different flux models was 
not significant. If we reject the AGN models having the parameters within the 3a contour, 
the critical values are fA,t ~ 0.4 and 6s t ~ 25° for the the 7-ray loud AGN models with 
UHECR flux proportional to the inverse square of the distance, and the critical values are 
fA.t ~ 0.7 and 6s,t ~ 20° for the the 7-ray loud AGN models with UHECR flux proportional 
to the 7-ray flux. This means that the models having large isotropic background with 7-ray 
loud AGN or the 7-ray loud AGN dominance models with large deflection by the intervening 
magnetic fields make them possible to describe the observed UHECR. At this stage, it is 
hard to confirm that 7-ray loud AGN are the source of UHECR. 
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